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On the interpretation of the apparent existence of a 
preferred magnetic polarity in extragalactic jet sources 
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ABSTRACT 

Contopoulos et al. recently argued that there is observational evidence for a pre- 
ferred sense of the Faraday rotation-measure gradients across jets from active galactic 
nuclei (AGNs). Such behaviour could arise if there were a deterministic relationship 
between the polarity of the poloidal magnetic field that threads the outflow and the 
sense of rotation of the outflow's source. Based on this interpretation, Countopoulos 
et al. suggested that their finding supports a model for the origin of cosmic mag- 
netic fields in a Poynting-Robertson process operating in AGN accretion discs. Here 
I point out that an alternative explanation of such a relationship could be that the 
Hall current plays a key role in the magnetohydrodynamics of the underlying disc. 
In this picture, the measured Faraday rotation is dominated by the contribution of a 
centrifugally driven wind that is launched from the weakly ionized outer region of the 
disc. Additional observations are, however, needed to verify the claimed behaviour. 
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1 INTRODUCTION 

In a recent paper, Contopoulos et al. (2009, hereafter 
CCKG09) compiled data on transverse gradients in the Fara- 
day rotation measure (FRM) across jets associated with 
active galactic nuclei (AGNs). The data were obtained by 
high-resolution (milliarcseconds, typically corresponding to 
a projected scale of parsecs at the source) multi-wavelength 
radio polarization observations. They found that the ma- 
jority (22/29) of sources in which gradients were detected 
relatively close to the base of the jet exhibited clockwise 
gradients. As explained by CCKG09, such gradients are ex- 
pected in outflows that contain a helical magnetic field (see 
Fig. 1 in Gabuzda 2008) . A helical field arises naturally when 
poloidal field lines that thread the outflow are twisted by the 
differential rotation of the source (an accretion disc or the 
central black hole in the case of AGNs), and the sense of the 
twist (as reflected in the FRM gradient) is directly related 
to the relative orientation of the poloidal field and the rota- 
tion vector (see Fig. 2 in Gabuzda et al. 2008). The inferred 
ordered magnetic field could originate in stars or interstel- 
lar gas and be dragged into the jet launching region by the 
accretion flow, or it could be generated in the disc by a dy- 
namo mechanism. In either case, one would expect an equal 
probability for the poloidal field component to have positive 
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or negative polarity. (For definiteness, I take the source to 
be a disc whose rotation vector points along the +z direc- 
tion of a cylindrical coordinate system {r, <j>, z} placed at the 
disc's centre, with the polarity of the magnetic field B de- 
fined as sgn{B z }.) CCKG09 suggested that the low (< 1%) 
probability for the compiled distribution to be a chance oc- 
currence indicates that there is a physical mechanism that 
acts to impose a positive magnetic polarity (i.e. B z parallel 
to the rotation vector) at the source of AGN jets. 

CCKG09 proposed that the relevant mechanism is the 
Poynting-Robertson 'Cosmic Battery' effect (Contopoulos 
& Kazanas 1998). In this picture, the radiation drag in- 
duced by the central AGN radiation field and acting pre- 
dominantly on the electrons in a rotating circumnuclear 
disc gives rise to an azimuthal current that generates a 
poloidal magnetic field with a unique (positive) polarity. 
Although the viability of this mechanism has been debated 
in the literature (e.g. Bisnovatyi-Kogan & Blinnikov 1977; 
Bisnovatyi-Kogan, Lovelace & Belinski 2002; Contopoulos, 
Kazanas & Christodoulou 2006; Christodoulou, Contopou- 
los & Kazanas 2008), CCKG09 argued that the indicated 
preponderance of positive polarities in AGN jets could be 
an important factor in its favour in view of the apparent 
difficulty of explaining this behaviour 'using any standard 
MHD model in the literature'. Here I point out that, in fact, 
a single- (positive) polarity disc outflow can naturally form 
also under 'standard MHD' conditions if the Hall current af- 



2 A. Konigl 



fects the magnetohydrodynamics (MHD) of the associated 
accretion flow. In this picture, the gas that dominates the 
observed FRMs corresponds to a wind that originates in the 
comparatively massive and weakly ionized outer region of 
the disc. If the claimed effect is real (which, in view of the 
considerable observational difficulties involved, requires fur- 
ther confirmation) and the disc-wind interpretation is cor- 
rect then the FRM measurements could yield valuable clues 
to the physical conditions in the outflow as well as in the 
underlying disc. 



2 HALL-CURRENT EFFECTS IN 
WIND-DRIVING DISCS 

For the purpose of illustration, I consider a simplified model 
of a wind-driving Newtonian disc that is in nearly Keplerian 
rotation (with azimuthal speed ~ vk) around a black 
hole of mass M, vertically isothermal (with sound speed 
c s ) and geometrically thin (so that the density scaleheight 
h(r) « r). I focus on systems that are axisymmetric and in 
a steady state on the dynamical time $1^ (r) = r/vn(r). I 
consider outflows that are launched in the form of a centrifu- 
gally driven wind (e.g. Blandford & Payne 1982), which are 
also relevant to jets that attain relativistic speeds (e.g. Vla- 
hakis &: Konigl 2004). The disc is assumed to be threaded 
by a large-scale, open magnetic field with an 'even' B z sym- 
metry about the mid-plane (corresponding to B r = — 
at z = 0). On going away from the mid- plane, the poloidal 
field lines bend away from the rotation axis, and the re- 
sulting radial component (B r ) is sheared by the disc's dif- 
ferential rotation to produce an azimuthal component (B^). 
The radial and azimuthal field components generated in this 
way satisfy BrB^ < 0. The B^ component induces a torque 
(oc rB z B$) that acts to brake the disc rotation. The mag- 
netic field removes angular momentum from the disc, and, if 
the field-line inclination at the disc's surface is large enough 
(B r /B z > 1/V3), a wind can be launched wherein the field 
transfers the angular momentum back to the gas (thereby 
accelerating it 'centrifugally'). 

The azimuthal shearing as well as the tendency of the 
accretion flow to advect the magnetic field lines inward are 
countered by the magnetic diffusivity of the gas. However, 
the diffusivity cannot be too large if vertical magnetic trans- 
port of angular momentum is to take place; in fact, the 
minimum-coupling condition is essentially the same as the 
corresponding condition for efficient radial magnetic trans- 
port through turbulence induced by the magnetorotational 
instability (MRI; see Konigl & Salmeron 2011). In the re- 
gion where this condition is satisfied, MRI may dominate at 
low disc elevations, with vertical transport by the large-scale 
field taking over further up (Salmeron et al. 2007, hereafter 
SKW07) . However, to simplify the discussion, I assume that 
the degree of coupling and the magnetic field strength are 
high enough for vertical transport to dominate already at the 
mid-plane. Wind-driving disc models of this type were pre- 
viously studied in the context of protostellar systems (e.g. 
Wardle & Konigl 1993; Konigl, Salmeron & Wardle 2010; 
Salmeron, Konigl & Wardle 2010; hereafter WK93, KSW10 
and SKW10, respectively). 

The vertical structure of the envisioned wind-driving 
disc region is shown schematically in Fig. 1. One can iden- 



tify three basic zones (WK93): the quasi- hydrostatic layer 
near the mid-plane (between z = and Zh = h) where the 
bulk of the matter is concentrated and most of the field-line 
bending takes place, a transition zone where the magnetic 
pressure comes to dominate the thermal pressure and the 
inflow gradually diminishes, and an outflow region (between 
the top of the disc at Zb and the 'sonic' critical surface at z B ) 
that corresponds to the base of the wind (with the mass out- 
flow rate effectively determined by imposing the regularity 
condition at z s ). The magnetic field lines are tied (possibly 
imperfectly) to the ionized disc component, and the colli- 
sional drag between this component and the neutrals trans- 
mits the magnetic torque to the bulk of the gas. Therefore 
fB0 < v$ inside the disc. The loss of angular momentum 
enables the neutrals to drift toward the centre, resulting in 
a radial drag on the ionized component that tends to impart 
an inward radial speed \vBr\ < \v r \ to the field lines. The ra- 
dial drag is balanced by the magnetic tension exerted by the 
outward-bent field lines, and this force, in turn, contributes 
to the radial support of the neutral gas inside the disc and 
causes it to rotate at sub- Keplerian speeds (v^ < vk)- Out- 
side the disc the situation is reversed: the field lines transfer 
angular momentum to the matter (vb4> > v$), which conse- 
quently rotates at super-Keplerian speeds (vj, > dk). 

The effect of the magnetic tension on the disc gas can 
be seen explicitly from the r component of the momentum 
equation, which in the quasi-hydrostatic region of a thin disc 
(in which the v z velocity component can be neglected and 
vertical gradients almost always dominate radial gradients) 
can be approximated by 
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(KSW10), where p is the mass density, c is the speed of 
light and the term involving the current density « 
(c/Air)dB r /dz represents the magnetic tension force. Under 
the same approximations, and assuming that angular mo- 
mentum transport by the large-scale field dominates, the <fi 
component of the momentum equation can be written as 
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where J r « —{c/A-K)dB^,/dz. Combining equations (1) 
and (2) yields 
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(where the subscript '0' denotes the mid-plane), which cor- 
roborates the inference that B r B^, < in the wind-driving 
region of the disc. 

The effect of the magnetic diffusivity on the disc struc- 
ture can be inferred from a consideration of Ohm's law, 
which for a weakly ionized gas takes the form 
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where 770, Vk an d Va are, respectively, the Ohm, Hall, and 
ambipolar diffusivities. The latter can be evaluated for an 
arbitrary charge composition, but it is often a good approx- 
imation to assume that the current carriers consist only of 
electrons (subscript 'e') and singly charged ions (subscript 
'i') with a mass ratio m c jra\ <g 1. In this case one can write 
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Figure 1. Schematic of the vertical structure of a wind-driving disc, showing a representative field line and the orientations of the 
Keplerian angular velocity fix and of components of the gas velocity (v), field-line velocity (vg), magnetic field (B), current density 
(J oc V X B) and magnetic force (oc J X B) vectors in the disc (2 < z^) and wind (2 > zj,) regions. See text for further details. 
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7/0 = c /4ira c , rjH = /3 e r]o and 7/a = /3 e ft?70 (Wardle 2007), 
where a c = n c e 2 /m c u c is the electron conductivity and 

eB 1_ 

rrijC Vj 

is the Hall parameter of a charged species j. In these expres- 
sions, n c is the electron number density (which is equal to 
rii), e is the electronic charge, = p < av >j n / (rrij + mn) is 
the momentum-exchange collision frequency with the neu- 
tral (subscript 'n') particles (where < av >j n is the charge- 
neutral momentum transfer rate) and B = |B| sgn{B z } is 
the signed magnetic field amplitude. The absolute value of 
Pj represents the ratio of the cyclotron frequency u C j to Uj 
and is a measure of whether the corresponding charged par- 
ticle is well (\/3j\ ^> 1) or poorly (\(3j\ -C 1) coupled to the 
magnetic field. Using rrii ~ 30 rrm (e.g. Draine, Roberge & 
Dalgarno 1983), the ratio of the ion and electron Hall pa- 
rameters can be estimated to be q = /?i//3 c ~ 1-3 x 10~ 4 \/T 
(where T is the gas temperature), which is typically <C 1. 

Note that, among the three diffusivity terms, only rju 
depends on an odd power of B - this is the origin of the 
dependence of the disc solutions on the magnetic field po- 
larity when the Hall current is dynamically important. To 
demonstrate the fact that negative-polarity field configura- 
tions are excluded in certain parameter regimes, consider 
the (f) component of equation (4) in the quasi-hydrostatic 
zone. As discussed in KSW10, ^ E^o = VB r oB z /c can 
be set equal to zero in this analysis without loss of gener- 
ality (see their Appendix A, where it is also argued that 
physical consistency requires one to consider all the terms 
in Ohm's law even in the nominally Hall-dominated regime). 
The advective term is approximately v r B z /c, where v r can 
be expressed in terms of J r using equation (2). Evaluating 
at z = 0, one gets 
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where the third expression is general and the fourth one 
applies in the case of a two-component plasma with q <C 1. 
In equation (6), the diffusivities are normalized using 77 = 
77//itCs, where /it = Cs/Qk is the tidal scaleheight, and the 
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the mid- plane ratio of the Alfven speed v\ = \B\/y/^Trp to 
the isothermal sound speed, and 
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the mid-plane ratio of the Keplerian rotation time to the 
neutral-ion momentum exchange time, are generally < 1 
and > 1, respectively (see KSW10). The two terms in the 
numerator of the third (or fourth) expression in equation (6) 
arise, respectively, from advective and Hall terms in Ohm's 
law. Their ratio is given by 
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where the last expression involves the Hall frequency 



1 

pn 



(9) 



(10) 



whose inverse represents the characteristic time-scale of the 
Hall effect in MHD (e.g. Pandey & Wardle 2008). (Here n is 
the particle number density, p the molecular weight and the 
subscript p denotes a proton.) When the ratio (9) is > 1 and 
sgn{B z } < 0, equation (6) implies that (dB r /dB^)o is > 0, 
which is inconsistent with the inference from the physical 
model (equation 3). Therefore, when this ratio is > 1, only 
positive-polarity disc outflows can form. 

Specializing to the case of a two-component plasma with 
1, the Hall-dominance condition can be written as 
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where the subscript 'd' denotes the disc, Xid = n^/na 
is the degree of ionization of the disc gas, ~ 2.33, 
M s = M/1O 8 M , S0.1 = Bo/0.1 G and r .i = r/O.lpc. 
When the parameter constraints on physically viable wind- 
driving molecular discs are analyzed in detail (see WK93, 
KSW10 and SKW10), one finds that the condition (11) is 
satisfied in two (out of the four identified) parameter regimes 
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in the nominal Hall diffusivity domain (|/3i| <C 1 <C |/3 e |) and 
in one regime (out of the three identified) in the nominal 
Ohm domain (|/3 e | <C 1). These are precisely the parameter 
regimes that admit only positive-polarity disc/wind solu- 
tions. (In the other viable-solution regimes both polarities 
are allowed.) In what follows, I consider whether the inequal- 
ity (11) can be satisfied in AGN discs but do not address 
the additional parameter constraints derived in the above- 
cited references, which could, however, be useful in a more 
detailed treatment. 



3 APPLICATION TO AGN DISCS 

Although the physical variables that appear in equation (11) 
cannot be directly determined from existing observations, 
one can nevertheless extract useful clues from the Faraday 
rotation measurements themselves. For example, Kharb et 
al. (2009) deduced from observations of parsec-scale jets in 
Fanaroff- Riley Type I (FR I) radio galaxies that the Fara- 
day rotation is likely associated with a helical (or toroidal) 
magnetic field and that it arises in a 'sheath' around the 
jet, which could correspond to a magnetized disc wind (e.g. 
Blandford & Levinson 1995). They inferred the magnetic 
field amplitude by equating it to the equipartition field 
(~ 5mG) derived from the synchrotron radio emission in 
the jet, which is reasonable if the jet is confined by the 
disc wind and both represent MHD outflows that are mag- 
netic pressure dominated. Using a representative FRM value 
of 200 rad mT 2 and a pathlength of ~ 1 pc through the 
sheath, they estimated a characteristic electron density of 
~ 0.03 cm -3 . This value could plausibly be ascribed to the 
atomic hydrogen layer that arises in the innermost zone of 
a photodissociation region formed when the AGN ionizing 
radiation (assumed not to be dominated by hard X-rays) im- 
pinges on the disc outflow (e.g. Tielens & Hollenbach 1985). 
A layer of this type is established just beyond the HII/HI 
transition zone (whose column density is < 10 20 cm -2 ) and 
has a characteristic electron fraction of ~ 3 x 10~ 4 (resulting 
primarily from singly ionized carbon atoms). This, in turn, 
implies a density n w « 10 2 cm" 3 (where the subscript 'w' 
denotes the wind) within the sheath. 

Now, if the parsec-scale sheath with \B\ < 10 mG can be 
identified with a centrifugally driven disc wind, the outflow 
likely originates on scales r < 0.1 pc, where the magnetic 
field strength is B(r) < 0.1 G, and its speed v w is > vk(t) 
(e.g. Blandford & Payne 1982). I adopt ro.i = 1 and Bo.i = 1 
as fiducial values and use Ms = 3 as a representative mass 
for radio- loud AGNs (e.g. Sikora, Stawarz & Lasota 2007). 
The mass outflow rate in the sheath can be estimated to 
be ~ r 2 /i w m p n w t> w k 8 x 10 23 gs _1 (using r w £s lpc, 
t) w ~ vk^'o.i = 1) and /i w ^ 1.27). Given that an efficient 
disc wind typically has a mass outflow rate that is about 
an order of magnitude smaller than the accretion rate (e.g. 
Konigl & Salmeron 2011), I adopt M 25 = M/10 25 gs" 1 « 1 
(corresponding to M ~ 0.16 Moyr -1 ) for the mass accre- 
tion rate at ro.i = 1. The mid-plane-to-disc-surface column 
density is given by N = M/4irrHdTn p \v r o\, and, by setting 
h « fiT (an upper limit due to the magnetic squeezing of the 
disc; see KSW10 and Fig 1), a lower limit on the mid-plane 
density can be obtained by writing n d0 = (2/n) 1 ' l2 N/hT- 
Using \v r o\ — 0.1eo.ic B (where eo.i is typically > 1; e.g. 



SKW07) and adopting T 500 = T /500K « 1 as a rough 
estimate of the mid-plane temperature at the radius of in- 
terest, I get N « 5 x 10 25 M 25 £o.i r o.l t mo' 2 cm ~ 2 and 
ndo > 3 x 10 11 (M 8 /3) 1/2 M 25 eo\ r^' 2 T M J cm" 3 . 

For the adopted values of Ms, ro.i and Bo.i, the condi- 
tion given by equation (11) can be fulfilled if Xido < 10~ 12 . 
The dominant ionization agent for the inferred values of 
r and N would likely be cosmic rays, and if one em- 
ploys the nominal Galactic values for the ionization rate 
(£cr ~ 1 x 10~ 17 s _1 ), the cosmic-ray attenuation grammage 
(x ~ 100gcrn~ 2 ) and the gas-to-dust mass ratio (w 100), 
one finds that the requirement on x id0 should be satisfied 
for the modelled disc (e.g. Umebayashi &: Nakano 1990). Al- 
though the charged-particle abundances depend on the un- 
certain heavy-element depletions and grain size distribution, 
ions and electrons have a significant influence on the mag- 
netic diffusivity under typical Galactic conditions so long 
as (h/^cr) exp (/idWpA^/x) is < 10 29 cm~ 3 s (e.g. Kunz & 
Mouschovias 2010). At higher values of this parameter com- 
bination, singly charged (+ and -) grains of equal mass be- 
come the dominant current carriers, which results in the 
Hall conductivity tending to zero (so Hall effects no longer 
play a significant role). The actual value of £cr in AGN 
discs is still unknown: on the one hand it might be substan- 
tially higher than in the Galaxy (e.g. Bayet et al. 2009), 
but on the other hand a super- Alfvenic disc outflow and the 
expected magnetic field geometry of the disc/wind system 
could strongly reduce the cosmic-ray flux reaching the mid- 
plane (e.g. Wardle 2007). In view of these uncertainties, the 
estimates employed in the evaluation of the inequality (11) 
need to be regarded as merely illustrative. These estimates 
nevertheless indicate that the proposed scenario is not im- 
plausible. Note that the wind- launching region is likely to 
contain magnetic fields of either polarity, both in the case in 
which the field is advected from a larger distance (for exam- 
ple, if it originates in stars with randomly oriented magnetic 
moments whose remnants are brought into the launching re- 
gion in sufficiently rapid succession) and in the case of a local 
disc dynamo (e.g. Blackman & Brandenburg 2003). There- 
fore, when the condition (11) is satisfied, a positive-polarity 
outflow should in general be able to form. 



4 DISCUSSION 

In the interpretation outlined in this Letter, the Faraday 
sheath that was inferred to be the likely source of the mea- 
sured parsec-scale FRM gradients is identified with a cen- 
trifugally driven wind that originates in a weakly ionized, 
molecular region of an AGN accretion disc at a distance 
> 0.1 pc from the central black hole. This picture is con- 
sistent with the idea that an MHD disc outflow constitutes 
the 'dusty torus' that underlies AGN unification schemes 
(Konigl & Kartje 1994). The suggestion that the gas in 
the outer regions of certain AGN discs is mostly molecular 
and dusty is supported by several observational findings. Al- 
though molecular line measurements do not yet resolve these 
regions, there are already strong indications that dense, ro- 
tating molecular gas is concentrated on scales > 10 pc 
around the centres of local AGNs (e.g. Hsieh et al. 2008; 
Hicks et al. 2009), and numerical simulations (e.g. Wada, 
Papadopoulos & Spaans 2009; Schartmann et al. 2010) indi- 
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cate that this gas likely forms a thin disc on smaller scales. 
Flattened distributions of dust have in fact been imaged in 
sources like NGC 1068 on > 1 pc scales by infrared obser- 
vations (e.g. Raban et al. 2009). Direct evidence for dense 
(n > 10 8 cm~ 3 ) molecular gas shielded from the central con- 
tinuum by a large dusty column (N > 10 23 cm' 2 ) has been 
provided in several AGNs (including NGC 1068; Greenhill & 
Gwinn 1997) by the detection of flattened distributions (im- 
aged on scales ~ 0.1 — 1 pc) of water masers that exhibit Ke- 
plerian (or nearly Keplerian) rotation-velocity profiles. The 
maser emission could arise either in the circumnuclear disc 
(e.g. Neufeld & Maloney 1995) or in dense clumps within 
a dusty disc wind (e.g. Kartje, Konigl & Elitzur 1999). 
Although all the above findings pertain to spiral galaxies 
(Seyferts and LINERs), there is also evidence that molec- 
ular gas is present in the centres of nearby radio-loud el- 
liptical galaxies (e.g. Ocafia Flaquer et al. 2010). Higher- 
resolution and more sensitive future observations, particu- 
larly with ALMA, should shed more light on this issue. 

A growing body of data indicates that relativistic AGN 
jets undergo the bulk of their acceleration and collimation 
on spatial scales not much smaller than ~ 0.1 pc (e.g. Sikora 
et al. 2005; Junor, Biretta & Livio 1999). Disc outflows could 
play an important role in these processes because efficient 
MHD acceleration to relativistic speeds requires lateral con- 
finement over extended spatial scales (e.g. Komissarov et al. 
2009) and because the presence of a surrounding disc wind 
could greatly facilitate the collimation of a relativistic jet, 
which might otherwise be hard to achieve (e.g. Gracia et al. 
2009) . It can be expected that the portion of the disc wind 
that is most effective in confining the jet on scales > 0.1 pc 
along the flow is launched from the disc surface at a compa- 
rable distance from the centre (see Spruit, Foglizzo & Stehle 
1997). This could naturally account for the proposed identi- 
fication of the Faraday-rotation sheaths of parsec-scale jets 
with magnetized winds that originate in the associated AGN 
discs at r > 0.1 pc. 
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